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ABSTRACT 

Context. Magnetic fields are at the heart of the observed stellar activity in late-type stars, and they are presumably 
generated by a dynamo mechanism at the interface layer (tachoclvne) between the radiative core and the base of the 
convective envelope. 

Aims. Since dynamo models are based on the interaction between differential rotation and convective motions, the 
introduction of rotation in the AT0N 2 . 3 stellar evolutionary code allows for explorations regarding a physically consis- 
tent treatment of magnetic effects in stellar structure and evolution, even though there are formidable mathematical 
and numerical challenges involved. 

Methods. As examples of such explorations, we present theoretical estimates for both the local convective turnover 
time (t c ), and global convective times (r g ) for rotating pre-main sequence solar-type stars, based on up-to-date input 
physics for stellar models. Our theoretical predictions are compared with the previous ones available in the literature. In 
addition, we investigate the dependence of the convective turnover time on convection regimes, the presence of rotation 
and atmospheric treatment. 

Results. Those estimates, as opposed to the use of empirically derived values of t c for such matters, can be used to 
calculate the Rossby number Ro, which is related to the magnetic activity strength in dynamo theories and, at least for 
main-sequence stars, shows an observational correlation with stellar activity. More important, they can also contribute 
for testing stellar models against observations. 

Conclusions. Our theoretical values of r c , r g and Ro qualitatively agree with those published by Kim & Demarque 
(1996). By increasing the convection efficiency, r g decreases for a given mass. FST models show still lower values. 
The presence of rotation shifts r g towards slightly higher values when compared with non-rotating models. The use of 
non-gray boundary conditions in the models yields values of r g smaller than in the gray approximation. 

Key words. Stars: evolution - Stars: interiors - Stars: rotation - Stars: pre-main sequence - Stars: magnetic activity - 
Stars: convection 



1. Introduction 

Magnetic activity in solar-type stars encompasses a variety 
of phenomena, such as star spots, activity cycles, heated 
outer atmospheres, X-ray emission, and many others. The 
driving mechanism for this activity is generally attributed 
to a dynamo that results from the interaction between ro- 
tation and convective motions in the star's outer envelope. 
Theoretical work by a number of researchers indicates that 
for main-sequence, solar-type stars the field is generated 
and amplified at the tachocline, the thin layer of differen- 
tial rotation between the convection zone and the nearly 
rigidly rotating radiative interior. For stars of a spectral 
type ranging from mid-F to early-M dwarfs, rotation and 
activity are thought to be controlled by this pro cess, also 
called an a—fl dynamo (|Mohantv fc Basril . I2003D . Its effi- 
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ciency is strongly dependent on the rotation rate and con- 
vective timescales. Young and rapidly rotating stars are, 
in general, very active. Specific models of dynamo theory, 
such as the a— 51 type, have been success ful in explaining 
the q ualitative features of solar activity (jWeiss fc Tobias! . 
I2000TI . 

Activity is strongly correlated with rotation velocity in 
the mid-F to mid-M dwarfs; it increases rapidly with the 
projected velocity, v sin i, then saturates above some thresh- 
old velocity (^10 km/s). This relationship is evident only 
down to K types; as we go from M to M6 types the rotation- 
activit y connection becomes less clear. iMohantv fc Basril 
(|2003l ) analyzed rotation velocities and chromospheric Ha 
activity, derived from high-resolution spectra, in a sample 
of mid-M to L field dwarfs. They found that, in the spectral 
type range M4-M8.5, the saturation-type rotation-activity 
relation is similar to that in earlier types, but the activ- 
ity saturates at a significantly higher velocity in the M5.5- 
M8.5 dwarfs than in the M4-M5 ones; this may result from 
a change in the dynamo behavior in later spectral types. 
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For fully convective stars, such as pre-main sequence 
late-type stars, this theory cannot be readily applied, as 
they miss the tachocline. However, since magnetic indi- 
cators such as active regions and strong flaring have also 
been reported for those stars, dynamo mechanisms operat- 
ing on the full convectio n region have also been proposed 
(e.g. iDurnev et afl . ll993h . 

On the other hand, observations of stellar activity in 
solar-type stars have shown a very tight relationship be- 
tween chromospheric Ca II H-K flux and the Rossby num- 
ber Ro = P m t/T c , where P rot is the rotation period and r c 
is the local convective turnover time. The Rossby number 
plays an important role in dynamo models, being related 
to the dynamo number D = Ro~ 2 which, in turn, is re- 
lated to the growth rate of the field. For solar-type stars, 
since t c cannot be directly measured, Ro is generally com- 
puted through a p olynomial fit of t c to the B— V color in- 
dex. For example, iNoves et alJ (|l984l ) give t c (B — V) as a 
theoretically-derived convective overturn time, calculated 
assuming a mixing length to scale height ratio a <~ 2. 

For pre-main sequence stars, Rossby numbers have 
been used to study the relations hip between the X-ray 
emission and magnet ic fields (e.g. iFlaccomio et al.l . l2003t 
iFeigelson et alJ . l2003f) . In this case, however, no clear re- 
lationship between activity and rotation is seen, as in the 
case of main-sequence stars, and one has to resort to evo- 
lutionary models for estimating Rossby numbers. 

Though current stellar evolutionary codes are not yet 
able to deal with magnetic fields, a first step towards this 
direction is the introduction of rotation on the models, as 
it is a key component of stellar dynamos. This is the case 
of the ATON 2 . 3 evolutionary code, in which both rotation 
and internal angular momentum redistribution have been 
introduced. Such capabilities allow us to make some ex- 
ploratory work towards a future version that can handle 
magnetic field generation from first principles. 

In this work we computed convective turnover times 
and Rossby numbers for a range of rotating low-mass stel- 
lar models, and discuss their behavior with time from 
the pre-main sequence to the zero-age main sequence. 
Our results are mainly compared with those found by 
iKim fc Demarquel (|1996l ). who provided the first self- 
consistent local and global convective turnover times 
and purely theor etical Rossby numbers (but see also 
I Jung fc Kim1[2007h . 

In Sect. [5] we describe the stellar evolutionary code used 
in this work (ATON code) as well as the specific inputs 
adopted for the present calculations. Results on Rossby 
number calculations are discussed in Sect. [31 The impact 
of other physical phenomena (like convection, rotation and 
atmospheric boundary conditions) is analyzed in Sect. 01 In 
Sect. [S]we present a simple application of our calculations 
to stars from the IC 2602 open cluster. Our conclusions are 
given in Sect. [5] 



temperatures; diffusive mixing and overshooting; and con- 
vection treatment under either the mixing length theory 
(MLT) or the Full Spectrum o f Tu rbulence (FST) from 
ICanuto fc Mazzitellil (jl99ll . ll992ft and lCanuto et alJ (jl996h . 
In addition to the standard gray atmosphere boundary con- 
ditions, the ATON code can also handle no n-gray atmo- 
spheri c integration based on the models by [ Heiter et al.l 
(|2002DHHauschildt et all (119991) : lAllard et all feoOOf) as de- 
scribed in Landin et al.l (|2006l ). 

The structural effects of rotati on were implemented in 
the ATON code according to the iKippenhahn fc Thomas! 
(|1970D method, as improved by lEndal fc Sofial (jl976fh 
Although the ATON code can handle the chemical mixing 
of species by using a diffusive approach (microscopic dif- 
fusion), the models presented here use instantaneous mix- 
ing, which is a good approximati on for pre-main sequence 
stars (jVentura fc Zeppieril . Il998al) . We acco unt also for the 
mixin g caused by rotational instabilities (jEndal fc Sofial 
I1978D . including dynamical instabilities (Solberg-H0iland 
and dynamical shear) and secular instabilities (secular 
shear, Goldreich-Schybert-Fricke and meridional circula- 
tion). Angular momentum redistribution in radiative re- 
gions was modeled through an advection-diffusi o n equ ation 
based on th e framework of IChabover fc Zahnl |l992) and 
IZahnl (|1992[ ). Angular momentum losses in the star's ex- 
ternal layers due to magnetized stellar winds are also taken 
into account in the form of a boundary con dition at the sur- 
face. W e adopted the prescription used in IChabover et aD 
(|1995f ) with a "wind index" n = 1.5, which reproduces well 
the lSkumamchl (|l972ft law v oc £~ 1/2 : 
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where w cr i t introduces a critical rotation level at which the 
angular momentum loss saturates (set to 5 Q@ in our mod- 
els, where fi© is the current solar surface rotation rate). 
The constant K in our models was calibrated by adjusting a 
1 M Q model so that its surface velocity matches the current 
solar rotation rate at the equator. M is the mass loss rate, 
which enters in Eqs. ([T]) and ^ in units of 10~ 14 M Q yr^ 1 , 
and which we set to 1.0. 

For the initial rotation rates, we adopt the relation be- 
tween initial angular momentum (J; n ) and stellar mass 
obtained from the correspondin g mass-r a dius and mass- 
moment of inertia relations from iKawalerl (|l987h : 



J in = 1.566 x 10 



5(1 
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Mo 
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gem s . 
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For more details about the treat ment given to rota t ion in 
the A TON code, see lMendesI (|l999h and lMendes et al.l (|l999l 
l200l . 



2. Input physics 

The ATON 2 . 3 code has many updated and modern fea- 
tures regarding the physics of stellar interi o rs, of w hich a 
full account can be found in IVentura et alJ (|l998bl ) . Some 
of its most important fea tures are: most up-to-date OPAL 



(llglesias fc Rogersl 19931) o paciti es, supplemented by those 



op 

of lAlexander fc Ferguson] (|1994l ) for lower (T < 6000iO 



3. Rossby number calculations 

Convective turnover times and Rossby numbers were com- 
puted for models ranging from 0.6 to 1.2 M (in 0.1 M© in- 
crements) with solar chemical co mposition. For ease ofcom- 
parison with a previous work by IKim fc Demarquel (|1996l ) , 
who provided theoretical calculations of Rossby numbers 
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for pre-main sequence stars, convection was treated accord- 
ing to the MLT, with the free parameter a set to 1.5 (which 
fits the solar radius at the solar age for a gray atmosphere, 
non-rotating model). Convection treatment and boundary 
conditions considerably affect the calculations of global con- 
vective turnover times and, consequently, the Rossby num- 
bers (see Sect. 14. ip . Rotation was modeled according to 
rigid body law in convective zones and loc al conservation 
of an gular momentum in radiative regions (|Mendes et all 

We followed the evolution of local and global convective 
turnover times (to be defined shortly) and Rossby numbers 
during the pre-main sequence for 0.6-1.2 M© stars and tab- 
ulated them together with the corresponding evolutionary 
tracks. Table[T] presents the 1 M Q models as an example of 
such tables. Column 1 gives the logarithm of stellar age (in 
years); Col. 2 the logarithm of stellar luminosity (in solar 
units); Col. 3 the logarithm of effective temperature (in K); 
Col. 4 the logarithm of effective gravity (in cgs); Col. 5 the 
logarithm of local convective turnover time (in seconds); 
Col. 6 the global convective turnover time (in days); and 
Col. 7 the Rossby number. 




10 6 10 7 10 8 10 9 10 10 10 1 

Age (yr) 

Fig. 1: Rotation period plotted against age, for each model 
mass. 
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Table 1: Pre-main sequence evolutionary tracks (including 
logr c , log T g and Ro) for 1 M© star a . 



log Age 
(years) 



log- 



6.2324 
6.4053 
6.5840 
6.7666 
6.9324 
7.0786 
7.2016 
7.2978 
7.3740 
7.4394 
7.5214 
7.9918 
9.0247 
9.3023 
9.4669 
9.5829 
9.6573 
9.6980 
9.7644 
9.8211 



log 1 off 

(K) 



logs 
(cgs) 



log T e 

(seconds) 



(days) 



Ro 



0.1468 
0.0217 
-0.1029 
-0.2172 
-0.2965 
-0.3278 
-0.2987 
-0.2231 
-0.1262 
-0.0387 
-0.0924 
-0.1590 
-0.1138 
-0.0862 
-0.0580 
-0.0288 
-0.0040 
0.0119 
0.0428 
0.0754 



3.6524 
3.6513 
3.6499 
3.6494 
3.6522 
3.6606 
3.6775 
3.6999 
3.7232 
3.7445 
3.7525 
3.7461 
3.7531 
3.7553 
3.7572 
3.7591 
3.7605 
3.7612 
3.7624 
3.7634 



3.855 
3.976 
3.095 
4.207 
4.297 
4.362 
4.401 
4.415 
4.411 
4.409 
4.494 
4.536 
4.518 
4.499 
4.479 
4.457 
4.438 
4.425 
4.399 
4.370 



6.3275 
7.1018 
7.0810 
7.0074 
6.9087 
6.7949 
6.6878 
6.5727 
6.4301 
6.2546 
6.1771 
6.1942 
6.1821 
6.1523 
6.1598 
6.1303 
6.1369 
6.1409 
6.1102 
6.1166 



766.6288 
488.2106 
369.1187 
291.4590 
227.1363 
179.1020 
134.8216 
99.6126 
72.7130 
53.4192 
41.6650 
46.0585 
42.4793 
41.0925 
39.4493 
38.1493 
38.0850 
37.5359 
36.8616 
36.2118 



0.1436 
0.0183 
0.0146 
0.0131 
0.0129 
0.0135 
0.0145 
0.0165 
0.0204 
0.0269 
0.0232 
0.0212 
0.2593 
0.7916 
1.0789 
1.4206 
1.5989 
1.7020 
2.0771 
2.3108 



"The complete version of the 
for all the masses of Table[2] 
form. 



table, including seven tracks 
will be available in electronic 



Figure [T] depicts the rotation period as a function of 
age for all models, and shows the typical spin-up during 
contraction followed by the longer phase of continuous spin- 
down. For the 1 M© model, this results in an initial velocity 
of nearly 3kms _1 at the beginning of the Hayashy phase; 
thi s is about one order of m agnitude below the value used 
bv iKim fc Demaraud (|l996h . When calculating convective 
overturn times, one must make an arbitrary assumption 
about where in the convection zone the dynamo is oper- 
ating. This assumption significantly affects the value of r c , 
since i t is strong l y dep endent on the depth. In this work, we 
follow [Gillilandl ( 1986f ) and calculate the "local" convective 
turnover time r c at a distance of one- half the mixing length, 
aHp/2, above the base of the convection zone. Its value is 
computed through the equation r c = aHp/v, where v is 
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Fig. 2: "Local" convective turnover time as a function of age 
for all models. 



the convective velocity. It should be noted, however, that 
there is no agreement in the literature reg arding the precise 
locati on where r c should be computed; iMontesinos et all 
(|200lD . for example, calculate it at distance of 0.95 Hp 
above the base of convective zone, with Hp taken as its 
value at the base of the convection zone. 

Figure^ shows the time evolution of the local convective 
turnover time (in seconds). For a given mass, it decreases 
during the Hayashy contraction and reaches its minimum 
when the contraction stops. After that, r c remains con- 
stant until the stars reach a main sequence configuration. 
However, most relevant for our purposes are Figs. [3]and[4j 
which show the profiles of the "global" convective turnover 
time and the "dynamo number", Ro~ 2 , respectively, as func- 
tions of age. The global convective turnover time (r g ) is 
defined as 



R h v ' 



(4) 
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Fig. 3: "Global" convective turnover time as a function of 
age for each model mass. 
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Fig. 4: The dynamo number as a function of age for each 
model mass. 



where Rb is the radius at the bottom of the convective zone, 
is the stellar radius and v is convective velocity. 
Figure [3] shows that r g follows the same behavior as the 
local convective turnover time, also decreasing substantially 
during contraction to the zero-age main sequence and, after 
that, remai ning nearly constan t and d epending only on the 
mass. As in lKim fc Demaraud (|l996f ). the local convective 
turnover time behaves like r g except for a scaling factor, 
because the convective turnover timescale is weighted to- 
wards the deepest part of the convection zone, where the 
shortcomings of the mixing length approximation are less 
important. With regard to Ro~ 2 , it is seen from Fig.0]that 
it follows r c during contraction but, after that, decreases as 
expected since the rotation period also increases. 

By using the evolutionary tracks, we constructed a set 
of isochrones for the ages of 0.2, 0.5, 0.7, 1.0, 2.0, 4.55 (solar 
age), 10 and 15 Gyr. In Table [2] we list their characteris- 
tics such as stellar mass in solar masses (Col. 1); logarithm 



Table 2: Isochrones for all models. 



Mass 






T s 




Rotation 


(M Q ) 


log (T«ff) 


log (L/L @ ) 


(days ) 


R - 2 


Period (d) 


0.2 Gyr 


0.60 


3.6148 


-1.13663 


121.1802 


8627.8177 


0.5581 


0.70 


3.6508 


-0.85884 


93.8181 


4711.2816 


0.5913 


0.80 


3.6878 


-0.60171 


77.0336 


4601.9900 


0.4777 


0.90 


3.7218 


-0.36391 


57.9367 


2093.6635 


0.5387 


1.00 


3.7480 


-0.14791 


43.8863 


1643.1720 


0.4534 


1.10 


3.7695 


0.05200 


30.5734 


875.9235 


0.4187 


1.20 


3.7908 


0.23895 


14.8697 


136.2882 


0.5509 


0.5 Gyr 


0.60 


3.6148 


-1.12952 


121.9241 


365.2809 


3.3443 


0.70 


3.6511 


-0.85206 


94.4531 


364.8090 


2.2956 


0.80 


3.6898 


-0.59006 


73.5144 


569.3372 


1.3836 


0.90 


3.7241 


-0.35056 


57.1289 


426.6247 


1.2034 


1.00 


3.7511 


-0.13146 


42.1241 


365.2252 


0.9081 


1.10 


3.7734 


0.07207 


28.4896 


268.6685 


0.7299 


1.20 


3.7938 


0.26043 


13.4697 


42.9852 


0.8630 


0.7 Gyr 


0.60 


3.6146 


-1.12785 


122.7990 


45.4688 


8.1315 


0.70 


3.6512 


-0.84948 


93.9007 


69.1059 


5.3534 


0.80 


3.6903 


-0.58622 


73.8996 


130.3651 


2.8620 


0.90 


3.7248 


-0.34541 


56.8609 


115.0898 


2.1514 


1.00 


3.7521 


-0.12448 


42.0428 


122.1356 


1.5846 


1.10 


3.7748 


0.08167 


26.7088 


92.5457 


1.1571 


1.20 


3.7950 


0.27275 


12.7590 


23.7916 


1.0832 


lGyr 


0.60 


3.6145 


-1.12602 


123.9869 


15.6947 


13.7087 


0.70 


3.6514 


-0.84667 


93.3632 


15.1359 


10.3630 


0.80 


3.6908 


-0.58179 


74.2237 


20.3958 


7.2070 


0.90 


3.7255 


-0.33908 


56.5470 


19.0868 


5.2591 


1.00 


3.7530 


-0.11547 


41.4909 


21.0411 


3.8355 


1.10 


3.7760 


0.09436 


25.8848 


15.9412 


2.8429 


1.20 


3.7962 


0.28993 


11.7507 


9.2329 


1.6338 


2 Gyr 


0.60 


3.6149 


-1.11995 


124.0207 


5.5953 


22.8118 


0.70 


3.6523 


-0.83796 


93.6465 


4.6988 


18.5225 


0.80 


3.6925 


-0.56838 


74.4090 


4.3392 


15.5263 


0.90 


3.7275 


-0.31932 


56.1457 


2.8531 


13.6315 


1.00 


3.7553 


-0.08640 


41.0063 


1.6047 


12.9745 


1.10 


3.7783 


0.13614 


24.9018 


0.6022 


12.7151 


1.20 


3.7979 


0.34593 


10.3207 


0.2400 


9.0570 


4.55 Gyr 


0.60 


3.6163 


-1.10689 


123.5031 


2.1587 


36.6458 


0.70 


3.6548 


-0.81720 


93.2665 


1.7031 


30.6947 


0.80 


3.6970 


-0.53436 


71.9753 


1.3278 


27.0554 


0.90 


3.7327 


-0.26633 


56.5844 


0.8790 


24.8868 


1.00 


3.7605 


-0.00376 


38.2314 


0.3903 


25.3930 


1.10 


3.7807 


0.25629 


21.7054 


0.1036 


27.8899 


10 Gyr 


0.60 


3.6196 


-1.08068 


121.1537 


0.8633 


55.7920 


0.70 


3.6610 


-0.77093 


90.9166 


0.6376 


47.6948 


0.80 


3.7070 


-0.45253 


69.9406 


0.4504 


43.8041 


0.90 


3.7428 


-0.12491 


51.0290 


0.2126 


44.6414 


15 Gyr 


0.60 


3.6228 


-1.05579 


119.3184 


0.5619 


69.1454 


0.70 


3.6674 


-0.72295 


89.8339 


0.3972 


60.0909 


0.80 


3.7163 


-0.35931 


66.2647 


0.2143 


58.7162 



of the effective temperature, in Kelvin (Col. 2); logarithm 
of the stellar luminosity, in cgs (Col. 3); global convective 
turnover time, in days (Col. 4); dynamo number (Col. 5); 
and the rotation period, in days (Col. 6). 

In Fig. [3] we show a plot of the global and local convec- 
tive turnover times versus log(T' ff ) for each age. Isochrones 
for the global convective turnover time versus period are 
shown in Fig. [51 The highest point of each curve corresponds 
to the lowest mass. For a given mass, r g varies slightly with 
the period. For a given period and age, the global convec- 
tive turnover time depends on stellar mass. 

Figure [7] shows the rotation period versus log(T c g). In 
that figure the rightmost points (those with the lowest 
temperatures) correspond to the lowest mass considered, 
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namely 0.6 M Q . We recall that in our rotating models con- 
vective regions are modeled by assuming rigid body rota- 
tion, while radiative regions rotate according to local con- 
servation of angular momentum. 

In Fig. [8] we plotted the inverse square of the Rossby 
number versus the rotation period. The lowest point of each 
line represents the highest mass. Figure shows the dy- 
namo number Ro~ 2 as a function of effective temperature 
and age. After establishing an empirical relation between 
Ro~ 2 and magnetic activity indices, one can use Fig. [S] to 
estimate the stellar mass and age from the effective tem- 
perature and an activity index. 

We believe that this set of results can be useful to sup- 
port observational studies of active pre-main sequence stars 
as well as for testing stellar models against observations. Of 
course, these determinations are model-dependent and sub- 
jected to the assumptions made in each model. 



W hen compared to the work of iKim fc Demaraiiel 
(|1996f ), our results concerning both the global and local 
convective turnover times agree well with theirs; but our 
Rossby number values are in general lower by one order 
of magnitude. Since Ro = P r ot/Tcj this difference can be 
mainly attributed to the different initial angular momen- 
tum of the models: while we use values of J\ u from E q. [3l 
which are different for each stellar mass. lKim fc Demarqud 
(|l996f ) adopted an initial rotation velocity of 30 km s _1 
for all their models, which results in higher initial rota- 
tion rates (for example, the initial rotation velocity for 
a 1 M Q model corresponds to near 3 km s _1 in our 
case). Some other factors can also contribute to the ob- 
served differences in Ro; for example, although convection 
is treated according to MLT approximation in both models, 
the convect ion efficiency is slightly dif ferent: while we adopt 
Q!mlt=1-5, IKim fc Demarquel (|1996f ) use q;mlt=1-86315; 
these values correspond to the a which reproduce the so- 
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lar radius at the solar age in e ach model. In addition , 
the opacities used in o u r mod els (jlglesias fc Rogersl fl99l 
lAlexander fc Ferguson , 1994) ar e mor e up-to-date than 
those used by Kim fc Demaraud (119961). n amely those by 
llglesias fc Rogersl (|199lD and lKuruczl (119911) . How ever, their 
atmospheric boundary conditions ( Kurucz , Il992f) are more 
realistic than the gray atmosphere approximation used in 
our work. Finally, the details concerning the internal redis- 
tribution of angular momentum can also be a source of such 
differences. 

Stellar activity is a complex phenomena rooted on the 
intricate interactions between magnetic fields, convection 
and differential rotation. In the case of the Sun and low- 
mass stars, it is generally accepted that magnetic fields are 
produced by a dynamo process inside the stars. 

Though stellar models th at include rotation , such a s the 
ATON 2 . 3 and the one from iKim fc Demarquel (|1996T) , can 
indeed compute purely theoretical Rossby numbers, it is 
clear that a more consistent scenario could be obtained if 
magnetic field generation was included in the models. For 
example, since the Rossby number Ro is closely related to 
the so-called "dynamo number" which plays a central role 
in the afi-dynamos, stellar models that include a dynamo 
process would allow a cross-check between these two values. 
This is however a difficult task, though some initiatives have 
been reported in the literature (e.g. iLvdon fc Sofial[l995l 
iLi et al.ll2006f ). As a matter of fact, work is in progress to 
include the effects of magnetic fields in the ATON 2 . 3 code 
itself, which we plan to address in a forthcoming paper. 

Asteroseismological data, which now are increasingly 
available through missions such as CoRoT and Kepler, 
are obviously of great importance to impose additional 
constraints on rotating stellar models, such as convective 
zon e depth, rotation rate and differential rotation (see 
e.g. iDupret et all 12004 IChristensen-Dalsgaardl 120081 and 
iTang et al.ll2008f) . 
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4. The impact of some physical phenomena on 
convective turnover times 

4.1. Convection 

In principle, a correct calculation of convective turnover 
times depends on a good knowledge about stellar convec- 
tion. Unfortunately, our poor understanding of this subject 
limits all determinations of such time scales. 

Presently, there are three main ways of computing 
convection in stellar e nvelopes: i) the tradi tional Mixing 
Length Theory (MLT, iBohm-Vitensel Il958f); ii) th e Full 
Spectrum of Turbulence (FST. ICanuto et al l Il996f ): and 
iii) MLT, in which the a value for each gravity and T e g 
is calibrated upon 2D or 3D hydrodynamical simulations 
(|Ludwig et al.l . ll999Ll2002l) . 

Since convective turnover times are strongly affected 
by the convection treatment, we analyze their dependence 
with different convection models and different convection 
efficiencies. With our version of ATON, we can compute con- 
vective turnover times by using two different convection 
regimes, the classical MLT and the FST. The other physi- 
cal parameters which do not concern convection, remain as 
described in Sect. [2] 

In Fig. [10] we show logr g as a function of stellar age 
for 0.5, 0.7, 1.0 and 1.2 M models computed with FST 
or MLT, with a=1.0, 1.5 and 2.0 for the latter case. After 



the radiative core is formed, the global convective turnover 
time decreases until the pre-main sequence contraction is 
over; the decrease of r g gets sharper as the stellar mass 
increases, reflecting the receding size of the convective re- 
gion. During the first stages of main sequence evolution r g 
remains nearly constant. For MLT models, the higher the 
convection efficiency (larger a parameter), the lower r g and 
hence a higher Rossby number; however, for masses above 
1 Mq, the situation is reversed as the star enters the main 
sequence. The FST models show still lower values of r g , as 
it describes the whole spectrum of convective eddies and 
consequently convective velocities are computed in an in- 
trinsically different way. 

The results obtained for the rotating, pre-main sequence 
models show the same behavior obta ined elsewhere for non- 
rotating, main sequence stars (e.g. iPizzolato et all l200lh 
regarding the chosen convection model. 

4.2. Rotation 



According to lDurnev fc Spruit] (|1979f ). rotation has a stabi- 
lizing effect on convection since it reduces the growth rate 
of instability. Even in the Sun, which has a low surface 
angular velocity, rotation can indeed influence convection. 

In Fig.[TTJwe show the global convective turnover time 
for 0.5, 0.7, 1.0 and 1.2 M Q , calculated with non-rotating 
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Fig. 10: Global convective turnover time as a function of 
age for selected models, for different masses and convective 
regimes. 



Fig. 1 1 : Global convective turnover time as a function of 
age for selected models, for different masses and rotation 
regimes. 



models (solid lines), with solid body rotation throughout 
the whole star (SB, dotted lines) and with a combination 
of solid body rotation in convective regions and differen- 
tial rotation in radiative zones (SB+Diff, dashed lines). 
Parameters related to other physical inputs remain as in 
Sect. [5] In the case of SB rotation, the whole star rotates 
with the same angular velocity during the evolution, while 
in the SB+Diff case our models take into account the sur- 
face angular momentum loss from stellar winds and the in- 
ternal redistribution of angular momentum. Since our mod- 
els correspond to slowly rotating stars (for example, our 1 
M models have an initial velocity of ~3 km/s at the begin- 
ning of the Hayashy phase) , the influence of rotation on the 
global convective turnover time becomes smaller for lower 
masses, being more evident at the end of the Hayashy phase 
and for stars with masses larger than 1 M Q . The values of 
r g obtained with models considering SB rotation are higher 
than those resulting from non-rotating models and models 
with SB+Diff rotation. Higher values of r g in the presence 
of rotation are indeed expected, since one of the main ef- 
fects of rotation is to mimic a lower mass star with a larger 
convective envelope. In addition, when comparing values 
of r g for SB rotation and for SB+Diff rotation, one would 
also expect lower values for the latter since in this case the 
angular velocity at the convection zone is lower due to the 
loss of angular momentum by winds and the redistribution 
of angular momentum. 

During the main-sequence evolution, the values of T g for 
SB+Diff rotation approaches the non-rotating ones as the 
rotation rate decreases, as can be seen in Fig. [IT] (dashed 
and solid curves). It is worth noting that these results 
consider the same initial angular momentum for SB and 
SB+Diff rotation, as given by Eq. ([3]). Some tests made 
with 1 M Q models reveal that r g increases with the initial 
angular velocity (and so the initial angular momentum). 

4.3. Non-gray atmosphere 

There is a good agreement in the literature that the use 
of non-gray atmosphere models as boundary conditions for 



stellar structure and evolution codes produces evolutionary 
tracks shifted toward lower effective temperatures in HR di- 
agram. The non-gray models include the treatment of atmo- 
spheric convection, which cannot be neglected at low T e ff's. 
The use of frequency-dependent opacities may also mod- 
ify the onset of convection within the atmosphere, so that 
the T ff of the evolutionary tracks can be strongly affected. 
Non-gray models produce stars with larger and deeper con- 
vective envelopes when compared with gray models, though 
their stellar radius is smaller. 

Here we report values of r g calculated with non-rotating 
models which use either gray or non-gray atmosphere 
boundary conditions. The other input physics are as in 
Sect. |31 In the first case the match between the interior and 
the external layers is made at the optical depth r=2/3, 
while in the non-gray treatment the bounda ry conditions 
are p rovided by atmosphere models from lAllard et al.l 
(|2000D , and a self-consistent integration is performed down 
to an optical depth t=10. 

The convective velocities v of non-gray models for the 
mass range considered show an interesting behavior when 
compared to the corresponding ones of gray models, they 
are lower during the pre-main sequence but are becoming 
higher as soon as they reach the zero-age main sequence 
(Fig. [T2"]) . This can be eas ily understood if one tak es in 
account that, as shown by lHarris et al.l (|2006l [2007) and 
now confirmed by our own models, the use of non-gray at- 
mospheres results in effectively larger opacities and lower 
effective temperatures during the pre-main sequence evolu- 
tion of low-mass stars. Sinc e a larger opac ity favors higher 
convective velocities (see e.g. lMaederl2009f) but a lower tem- 
perature does the opposite, our models shows that this lat- 
ter effect prevails over the former during the pre-main se- 
quence. From the zero-age main sequence on the situation 
reverses, as the temperatures become nearly the same for 
both non-gray and gray atmosphere models, but the opac- 
ity remains larger for the non-gray models, and so results 
in higher convective velocities. 

These general effects on the convective turnover time 
can be seen in Fig. [121 where T g is plotted against the stel- 
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Fig. 13: Global convective turnover time as a function of 
age for 0.5, 0.7 1.0 an 1.2 M models (from top to bottom) 
and for different atmospheric treatment. 



lar age for 0.5, 0.7, 1.0 and 1.2 M Q models computed with 
gray (solid lines) and non-gray (dotted lines) atmosphere 
boundary conditions. Curves at the top of the figure cor- 
respond to the 0.5 Mq, and the global convective turnover 
time decreases as the stellar mass increases. 



5. Applications and comparisons with observations 

Star cluster data have particular importance in testing evo- 
lutionary models with rotation. The spin-down of young 
cluster stars places a lower limit to the time scale for the 



angular momentum redistribution at the base of the sur- 
face convection zone. Transport of the angular momentum 
leads to rotationally induced mixing which reproduces the 
observed Li depletion in the Sun an d solar analogs in open 
clusters (|Pinsonneault et"al~l . I1989T ) ■ Observations of rota- 
tion rates of stars in clusters of different ages help us to draw 
a scenario for the stellar angular momentum evolution. The 
broad period distribution of zero-age main-sequence stars 
and the existence of slow rotators can only be explained 
if during early phases of their evolution some stars lose 
a significant fraction of their angular mom entum which 
in ad dition must be different for each star (|Lamm et al.l 
l2005f ) . A number of studies have been published in the lit- 
erature showing how observational data from open clusters 
can impose constra ints on stellar models as for e x ample 
lLandin et al.l (|2006l ) and iRodrieuez-Ledesma et all (|2009f) 
for the Orion nebula cloud or lMeibom et al.l (|2009t) for M35. 

In this context, and to give a simple application to our 
results, we used the AT0N code to calculate Rossby numbers 
for a representative sample of solar-type stars in the young 
(30 Myr) open cluster IC 2602, located at the southern 
hemisphere at a distance of about 150 pc. 

As can be recalled from the introduction, the Rossby 
number is not a quantity directly obtained from observa- 
tions, since it is the ratio of the local convective turnover 
time t c to the rotation period P ro t- While r c can be de- 
rived fro m evolutiona r y mo dels or through a polynomial fit 
to B-V (|Noves et allfl98l . P ro t can be obtained observa- 
tionally or computed through rotating evolutionary mod- 
els, though this last method is obviously model-dependent. 
In this section, we computed Ro for IC 2602 stars in two 
ways: (i) by using r c from our models and P rot from observa- 
tions and (ii) by using both r c and P rot calculated through 
our models. We designate the Ro calculated as described 
in case (i) as "semi-theoretical" Rossby numbers and those 
computed as in case (ii) as "purely theoretical" Rossby num- 
bers. 

The semi-theoretical values of Ro were obtained by 
using the ob s erved rotation periods of IC 2602 stars 
(|Barnes et all Il999h . which are given in the range from 
0.2 days to ~10 days as a function of the B— V color index, 
and t c calculated by our models at the age of IC 2602. Since 
our models provide the local convective turnover time as a 
function of effectiv e temperature and gravity, we used the 
iBessel et al.l (|1998l ) relations, which provide B— V for a grid 



of T c fj and log g, to obtain r c as a function of B— V. In this 
way, by using t c {B — V) from the models and P rot (P — V) 
from observations, we calculated Ro(B — V). The purely 
theoretical Rossby numbers, in turn, were obtained from 
rotating models (SB+Diff) with gray atmosphere boundary 
conditions, MLT convection treatment (a=1.5) and inter- 
nal angular momentum redistribution with surface angular 
momentum loss, for the mass range of 0.2-1.4 M at the 
age of IC 2602. Those models provide us with r c and P ro t 
as a functio n of T e g and gravity, and, again by means of 
IBessel et al.l fl998) relations, we obtained the purely theo- 
retical Rossby numbers corresponding to the observed B— V 
of IC 2602 stars. In Fig. [T3] we present our purely theo- 
retical Rossby numbers (circles) and semi-theoretical ones 
(triangles) for the IC 2602 star sample. Cros ses in Fig. [T4"l 
repres ent the Rossby numbers presented by iBarnes et al.l 
(l!999h and were plotted here for comparison purposes. 
Thes e Ro values were c ompu ted with the use of models 
from iKim fc Demarauel (|l996h and are also designated as 
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"semi-theoretical", since their way of computing Ro is sim- 
ilar to ours. 

For a better comprehension of differences in these Ro 
calculations, we also display the corresponding values in 
Tabled in which Col. 1 gives the identification of each star 
of the IC 2602 sample; Col. 2 the color index B-V; Col. 
3 the purely theoretical Rossby number; Col. 4 the semi- 
theoretical Rossby number; a nd Col. 5 the s emi-t heoretical 
Rossby number presented by iBarnes et al.l (|1999f ) . 

When com paring our semi-the oretical Rossby numbers 
to those from IBarnes et all (|1999t ). the differences are ob- 
viously due only to the corresponding local convective 
turnover times, since the rotation periods are the same. 
As discussed in Sect. [31 these discrepancies in r c can be 
attributed to differences in the convection efficiency, at- 
mospheric boundary conditions, opacities and treatment 
of transport of angular momentum. Besides, although the 
distributio ns of our sem i -theor etical Rossby numbers and 
those from I Barnes et aD (|l999h are quite similar in Fig. [T4l 
despite the differences of r c presented in Tabled one 
must be cau t ious when interpreting those distributions. 
IBarnes et al.l (|l999f ) only mentioned that they used the r c 
values from the models of iKim &; Demarqud (|1996l ) to cal- 
culate Ro for the individual stars of their sample, but they 
do not disclose how they obtained their values of r c as a 
function of B— V. 

Our purely theoretical Ro values, in turn, behave very 
differently from the semi-theoretical ones due to the differ- 
ent way in which rotation periods are obtained in these two 
cases. In the first case the rotation periods start from an 
initial value of the angular momentum given by Eq. and 
evolve according to local conservation of angular momen- 
tum in radiative zones and rigid body rotation in convective 
ones; this initial value depends on the stellar mass, i.e. also 
on B— V. On the other hand in IC 2602, as in most young 
open clusters, there exist simultaneously ultrafast, interme- 
diate and slow rotators independently of the stellar mass; in 
order to cover this broader period distribution we need to 
consider some disk regulation mechanism to describe the 
magnetic coupling of the central star to its circumstellar 
disk. The role played by this disk regulation in the rota- 
tional evolution of young clusters and its effects on Rossby 
numbers will be discussed in a forthcoming paper. 

6. Conclusions 

Our results show the same trends for the Rossby number Ro 
an d the local convective turn over time r c found in the work 
by IKim fc Demarqud (|l996f ) as, for example, the decrease 
of t c during the pre-main sequence phase and its nearly 
constant value from that point on. Despite this, our values 
of Ro are in general lower by one order of magnitude, mainly 
due to differences in the initial rotation rates. 

As already expected, r g , the global convective time, de- 
creases as convection efficiency increases, and consequently 
Ro increases; the opposite situation occurs when stars with 
masses larger than 1 M Q enter the main sequence. We also 
found that the values of r g obtained for FST models are 
lower than those obtained with MLT (a— 1.0, 1.5 and 2.2), 
since the former compute the convective velocities in an 
intrinsically different way, which describes the whole spec- 
trum of convective eddies. 

The effect of rotation on r g is less important for masses 
lower than 1 M Q , but for masses larger than this threshold 
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Table 3: Rossby numbers and B— V for IC 2602 stars. 



Name 3 


B-V a 


Ro 


Ro 


i?o a 




(Pur. Theor.) 


(Semi- Theor.) 


B134 


0.95 


0.0132 


0.1767 


0.110 


W79 


0.83 


0.0176 


0.2107 


0.117 


R15 


0.93 


0.0136 


0.0967 


0.059 


R24A 


1.43 


0.0039 


0.0127 


0.013 


R27 


1.50 


0.0023 


0.0328 


0.045 


R29 


1.11 


0.0097 


0.0440 


0.031 


R43 


0.95 


0.0132 


0.0203 


0.013 


R52 


1.07 


0.0105 


0.0084 


0.006 


R56 


1.43 


0.0039 


0.0418 


0.043 


R58 


0.65 


0.0281 


0.0355 


0.015 


R66 


0.68 


0.0261 


0.1748 


0.082 


R70 


0.69 


0.0255 


0.2212 


0.107 


R72 


0.64 


0.0288 


0.0691 


0.029 


R77 


1.47 


0.0030 


0.0855 


0.104 


R83 


0.62 


0.0303 


0.1226 


0.049 


R88A 


1.20 


0.0080 


0.0035 


0.003 


R89 


1.24 


0.0073 


0.0757 


0.059 


R92 


0.67 


0.0267 


0.1117 


0.050 


R93 


1.37 


0.0050 


0.0802 


0.075 


R94 


1.39 


0.0047 


0.0295 


0.028 


R95A 


0.87 


0.0154 


0.0362 


0.021 


R96 


1.25 


0.0071 


0.0282 


0.022 


a [Barnes 


et al.l (1999) 









the influence of rotation on r g is more evident. Solid body 
(SB) rotation produces stonger deviations on r g relative to 
the non-rotating value than the SB+Diff case. By increasing 
the initial angular momentum, r g tends to increase as well. 

Among the effects analyzed here, the ones yielded by 
atmosphere boundary conditions are those that have the 
most influence on the values of r g . Models which use non- 
gray boundary conditions produce values of r g which are 
lower than those of the gray ones. 
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Our models were applied to calculate Rossby numbers 
as a function of the B— V color index for a sample of stars 
from the IC 2602 open cluster. Semi-theoretical values of Ro 
calculated wi th our r r . are , on av erage, higher than those 
presented by iBarnes et al.l (|1999f ) . Purely theoretical and 
semi-theoretical Rossby numbers have a different behav- 
ior due to the different origin of the rotation period which 
composes Ro in each case. 
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